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In contrast to gravity in the weak-field regime, which has been subject to numerous experimental 
tests, gravity in the strong-field regime is largely unconstrained by observations. We show that 
gravity theories that cannot be rejected by solar system tests but that diverge from general rel- 
ativity in the strong-field regime predict neutron stars with significantly different properties than 
their general relativistic counterparts. In particular, the range of redshifts of surface atomic lines 
predicted by such gravity theories is significantly larger than the uncertainty introduced by our lack 
of knowledge of the equation of state of ultra-dense matter. Measurements of such redshifted lines 
with current X-ray observatories such as Chandra and XMM-Newton can thus provide interesting 
new constraints on strong-field gravity. 



The rapid variability and high luminosities of accreting 
compact objects strongly suggest that their high-energy 
emission originates very deep in their gravitational poten- 
tials. As a result, the properties of their X-ray and 7-ray 
emission can, in principle, be used to test directly the 
strong-field regime of a gravity theory. Unfortunately, 
there are many astrophysical considerations that com- 
plicate the modeling of the emission of accreting neu- 
tron stars and black holes and it is difficult to disentan- 
gle the effects of gravity from those due to astrophysi- 
cal processes. For this reason tests of gravity involving 
high-energy phenomena around compact objects have re- 
ceived little attention in the past. This is despite the fact 
that the behavior of gravity in the strong-field regime is 
largely unconstrained by observations. 

The properties of a neutron star and its external space- 
time, in contrast to the high-energy phenomena above its 
surface, are determined solely by the strong-field behav- 
ior of gravity and the equation of state. To date, studies 
of neutron stars have focused on better determining the 
equation of state, whose details are still debated. We will 
argue, however, that the uncertainty in the behavior of 
strong-field gravity introduces much greater differences 
in neutron-star properties than do current uncertainties 
in the neutron-star equation of state. Furthermore, new 
observational developments allow for the measurement of 
the gravitational redshift of surface atomic lines, which 
is the cleanest source of information about the mass, ra- 
dius, and spin of neutron stars. The first such obser- 
vations, with Chandra and XMM-Newton, have already 
been reported U 01 ■ 

Testing rigorously General Relativity (GR) in the 
strong-field regime requires a general framework of 
which GR is a subset; this would be equivalent to the 
Parametrized Post-Newtonian formulation in the weak- 
field. Such a framework is yet to be constructed. Thus, 
in order to provide a proof of principle for the efficacy 



of the proposed tests, we will investigate here the con- 
straints that can be imposed on a parametrized subclass 
of scalar-tensor theories (see, e.g., Ref. Q)- As an ad- 
ditional example, we will also consider the limits that 
can be placed on Rosen's bimetric theory although 
binary-pulsar experiments have already excluded the bi- 
metric theory [5j, our results demonstrate the potential 
of the proposed tests for constraining gravity theories. 

We consider first scalar-tensor theories. The general 
class of scalar-tensor theories, where the gravitational 
force felt by matter is mediated by both a rank-two ten- 
sor, g^^, and a scalar field, 0, is one of the most natural 
extensions of Einstein's theory. The action is Q 
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Here V&m refers collectively to all matter fields other than 
(j), G* is a dimensional constant, and A((j)) is a function 
of (j). The tensor g^^i, is the metric in the Einstein frame 
and the Brans-Dicke frame tensor is g^^, — A'^{ip)g^. 
For A{(p) unity, we recover GR. 

Solar system experiments constrain the coupling func- 
tion A to be very fiat at the cosmological value of the 
scalar field cj). If, however, the scalar field inside a com- 
pact object fiuctuates far enough away from its cosmo- 
logical value to discover a steep part of A, the system 
will be able to reach the more energetically favorable 
configuration with large non-zero (j) near the center of 
the object. This non-perturbative effect leads to neutron 
stars that are significantly more massive and larger than 
their GR counterparts, but it is invisible to solar sys- 
tem experiments that probe the weak-field regime. This 
is the phenomenon of "spontaneous scalarization" dis- 
covered by Damour & Esposito-Farese G], who used the 
Hulse- Taylor pulsar to put constraints on the strong-field 
regime of such theories. 
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Following Damour & Esposito-Farese |g, we choose 
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where /3 is a real number. Spontaneous scalarization oc- 
curs only when (3 < —4.85. Throughout this paper we 
will be interested only in solutions where 0o, the cosmo- 
logical value of the scalar field, is zero and v4(0o) = 1- 

Following Ref. 0|, we write the relativistic equations 
of stellar structure for a non-rotating neutron star in 
scalar-tensor gravity. We write the metric as 
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and describe the matter fields as a perfect fluid in the 
physical frame, i.e.. 



Tf^f = ip + p)up.Ui, + pg^^u- 



(4) 



Here p is the pressure, p is the density, and is the 
four-velocity of the fluid. The differential equations to 
be integrated are then 



v' = 87rg, ^ + r^^ + ^ 
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in Cook et aL[9|. In order of increasing stiffness, these 
are EOS A [ig and EOS UU [ill. Neutron star mod- 
els computed with these two equations of state bracket 
the uncertainty introduced by our inability to calculate 
from first principles the properties of ultra-dense matter, 
when condensates or unconfined u-d-s quark matter is 
not taken into account. 

Having chosen an equation of state, we then integrate 
the system of equations (|K|- l(Tn|l from the center of the 
star, where we specify p(W) = i/(0) = A^(0) = 0, 0(0) — 
4>c, and p(0) = pc, to the surface, where p{R) = 0. We 
then integrate equations from the surface of the 

star to infinity to determine the form of the metric exte- 
rior to the star. We note that the Arnowitt-Deser-Misner 
(ADM) mass felt by an observer far away is in general 
different to the mass that contributes to the redshift. 

In addition to considering the scalar-tensor theory, we 
also investigate the bimetric theory, proposed by Rosen 
0. This involves, in addition to the dynamical metric 
f/^i/, a nondynamical flat metric rj^^. 

The metric g'^i/ in Rosen's theory may be written as 



ds' 



,2 /) I .2 



% (11) 



where m$(r) and m!^(r) are given, interior to the star, 
by 



m^{r)='i7r e'^+^^ip + Sp^dr, (12) 

^0 



TOA(r)=47r/ e'^+^{p-p)r^dr, (13) 
Jo 

and $ and A are given by the field equations 

$' = GTO<i.(r)/r2, (14) 
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Here N is the baryon number, n is the number density, 
and primes denote derivatives with respect to r. We sup- 
plement these equations with a (zero temperature) equa- 
tion of state, p = p{p) and n ~ n{p). 

We choose two commonly used equations of state, 
which cover a broad subset of the wide range discussed 



A' = -GTOA(r)/r2. 
The equation of hydrostatic equilibrium is 
p' = -(p + p)^'. 
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Outside the star, m^{r) = m^{R) = M$ and m{^{r) = 
mtv{R) = Ma are constants. M$ is equal to the Kepler- 
measured mass at large distances. We specify an equa- 
tion of state, p — p{p), and the boundary conditions 
/j(0) = /5c, $((»)= A(oo) = 0. 

Figure ^ shows the relationship between ADM mass 
and neutron-star radius for the three theories under con- 
sideration. For the case of the scalar-tensor theory, we 
show, as an example, the case [3 — ~8, which is com- 
parable to the most negative value of /3 not yet ruled 
out by binary pulsar timing For each of the three 
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FIG. 1: Mass-radius relations for neutron stars in the GR, bi- 
metric, and scalar-tensor theories. For the scalar-tensor the- 
ory, we plot the relation for only one value of the parameter 
/3. The shaded regions show the allowed values of Af and 
R for equations of state with stiffness between EOS A and 
EOS UU. 

theories, we plot a shaded region corresponding to the 
allowed values of mass and radius for equations of state 
with stiffness between EOS A and EOS UU. We include 
only the stars stable to radial perturbation and, for con- 
figurations in the scalar-tensor theory we show only the 
one energetically preferred. 

As can be seen immediately, all three theories allow 
neutron stars with masses in the currently observed range 
of 1.35 — I.SMq to exist a result, measurement 

of neutron star masses alone cannot distinguish between 
them. This degeneracy is broken, however, by consider- 
ing the predicted radii, fn the astrophysically relevant 
range Madm > I.BA/q, the three theories occupy mutu- 
ally exclusive regions in the (AfADM,^) space. Observa- 
tions that can put limits on both the mass and radius 
of a neutron star can thus constrain, without ambiguity, 
the permissible set of gravity theories. 

One such observation is of the surface redshift fac- 
tor, z, relative to infinity. This is defined by E^c = 
[1/(1 + z)]Esurf, where E'surf is the energy of a photon 
emitted from the stellar surface and measured at infin- 
ity with energy i?oo- In the GR Schwartzschild met- 
ric, this provides a direct measurement of only the ratio 
Madm/R- In the case of the scalar-tensor and bimet- 
ric theories the relationship between Madm, R, and z is 
more complicated. In all cases, however, the formula, 

E^=(9^Y\^^,,, (17) 
V 500,00/ 

holds, where 500, s is evaluated at the surface of the star 
and 500,00 is evaluated at infinity. 

Figures[21and|21show the redshift factor z for the scalar- 
tensor and GR theories as contours on a plot of the mass 



FIG. 2: Contours of constant redshift in the scalar-tensor 
theory, as a function of Madm and /3, for EOS A. The dashed 
line shows the value of Madm, as a function of (3, for which 
the baryonic mass is equal to IAMq. The heavy lines on the 
plot show the mass limits for neutron stars and the bounding 
region where spontaneously scalarized stars are produced. No 
stable neutron stars in the crosshatched region exist. 
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FIG. 3: Contours of constant redshift in the scalar-tensor 
theory, as a function of A/adm and (3, for EOS UU. See figure|21 
for details. 

Madm versus the parameter /?, for EOS A and EOS UU, 
respectively. As discussed above, the phenomenon of 
spontaneous scalarization occurs only for /3 ^ —4.35, and 
so above this value the contours are parallel to the (3 
axis and equal to the GR value. The heavy lines on the 
plot show the upper mass limits for neutron stars as well 
as the bounding region, in the (Madm,/?) space, where 
spontaneously scalarized stars are produced. The dashed 
line shows the value of Madm as a function of P for which 
the baryonic mass is equal to 1.4M©. Note the z — 0.35 
contour, the value recently measured by Cottam et al. 
in the neutron star source EXO 0748-676. 
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FIG. 4: Mass-redshift relations for the GR and bimetric the- 
ories. The dashed line shows z = 0.35, the value reported by 
Cottam et al. 

The redshifts found for stars in the scalar-tensor the- 
ories differ greatly from their general relativistic values 
for both equations of state. The trend is for the redshift 
values to increase (nearly always) monotonically as /3 be- 
comes more negative. If we can place constraints on the 
neutron star mass, a redshift measurement will put an 
upper bound on the value of — /3. 

Figure 01 shows the redshift factor z as a function of 
A^ADM for neutron stars produced in the GR and bimet- 
ric theories. We plot two shaded regions, one for GR 
and one for the bimetric theory, that cover the range be- 
tween EOS A and EOS UU. The dashed hue shows the 
recent z — 0.35 measurement 0. Even when a GR and 
a bimetric star have the same redshift, the mass serves 
to break the degeneracy, providing a test of strong field 
gravity despite the uncertainty in the equation of state. 

The recent measured redshift of z = 0.35 \\ cannot 
by itself constrain strong-field gravity in the theories we 
consider here. However, strong limits can be placed if 
the mass of the neutron star is measured directly from 
binary dynamics or constrained by arguments from for- 
mation mechanisms (e.g., restricting the baryonic mass 
to be above the Chandrasekhar limit for a degenerate 
core.) The neutron star EXO 0748-676, for which this 
redshift was measured, is a member of an eclipsing binary 
system with an orbital period of 3.8 h [3| and hence is 
a prime candidate for producing a mass measurement 
that can provide such a constraint. Simply requiring the 
baryonic mass of the star to be greater than 1.4Mq (the 
Chandrasekhar limit for its degenerate progenitor) places 
an upper limit on -/3 of 7 (EOS A) to 9 (EOS UU). As 
expected, the constraint on strong-field gravity depends 
only weakly on the equation of state. 

In addition to the redshift of atomic lines, a number of 
other observational characteristics of neutron stars de- 
pend - in more complicated ways - on the of nature 



strong-field gravity. These are the maximum allowed 
spin frequency, the Eddington limit for accreting stars, 
bursting behavior, cooling rates, and the frequencies of 
quasi-periodic oscillations. The constraints imposed by 
quasi-periodic oscillations will be reported in a forthcom- 
ing paper. 

The constraints on scalar-tensor discussed here have a 
wider applicability than the structure of compact objects 
because varieties of scalar-tensor theories naturally arise 
as low-ener gy l imits of gravitational theories in higher 
dimensions |l7| . They are also under consideration in 
cosmology as explanations of evidence for an accelerat- 
ing universe. To date, little attention has been given 
to the implications of such theories for the properties of 
compact objects that are now observable. The uncertain- 
ties in the equation of state at ultra-high densities were 
thought to make such studies fruitless. Our investigation 
has shown that this is not the case and that measure- 
ments of the redshifts of surface atomic lines from neu- 
tron stars can provide new tests of strong-field gravity. 

We thank A. Loeb, J. Maldacena, C. Miller, and F. 
Ozel for useful discussions. SD acknowledges the support 
of an NSF Graduate Research Fellowship. DP acknowl- 
edges the support of NSF grant PHY-0070928. 



[1] J. Cottam, F. Paerels, and M. Mendez, Nature 420, 51 
(2002). 

[2] D. Sanwal, G. G. Pavlov, V. E. Zavlin, and M. A. Teter, 

Astrophys. J. 574, L61 (2002). 
[3] T. Damour and G. Esposito-Farese, Class. Quantum 

Grav. 9, 2093 (1992). 
[4] N. Rosen, J. Gen. Rel. and Grav. 4, 435 (1973). 
[5] C. M. WiU and D. M. Eardley, Astrophys. J. 212, L91 

(1977). 

[6] T. Damour and G. Esposito-Farese, Phys. Rev. Lett. 70, 
2220 (1993). 

[7] C. Will, Theory and Experiment in Gravitational Physics 

(Cambridge University Press, 1993). 
[8] T. Harada, Phys. Rev. D 57, 4802 (1998). 
[9] G. B. Cook, S. L. Shapiro, and S. A. Teukolsky, Astro- 
phys. J. 424, 823 (1994). 
[10] V. R. Pandharipande, Nucl. Phys. A 174, 641 (1971). 
[11] R. B. Wiringa, V. Fiks, and A. Fabrocini, Phys. Rev. C 

38, 1010 (1988). 
[12] T. Damour and G. Esposito-Farese, Phys. Rev. D 54, 
1474 (1996). 

[13] S. E. Thorsett and D. Chakrabarty, Astrophys. J. 512, 
288 (1999). 

[14] O. Barziv, L. Kaper, M. H. Van Kerkwijk, J. H. Telting, 

and J. Van Paradijs, a 377, 925 (2001). 
[15] J. A. Orosz and E. Kuulkers, MNRAS 305, 132 (1999). 
[16] A. N. Parmar, N. E. White, P. Giommi, and 

M. Gottwald, Astrophys. J. 308, 199 (1986). 
[17] J. M. Overduin and P. S. Wesson, Phys. Rep. 283, 303 

(1997). 



